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ABSTRACT 

The chemical strategies used by ribozymes to 
enhance reaction rates are revealed in part from 
their metal ion and pH requirements. We find that 
kinase ribozyme K28(1-77)C, in contrast with previ- 
ously characterized kinase ribozymes, requires Cu^^ 
for optimal catalysis of thiophosphoryl transfer from 
GTPyS. Phosphoryl transfer from GTP is greatly 
reduced in the absence of Cu^^, indicating a 
specific catalytic role independent of any potential 
interactions with the GTPyS thiophosphoryl group. 
In-line probing and ATPyS competition both argue 
against direct Cu^^ binding by RNA; rather, these 
data establish that Cu^^ enters the active site 
within a Cu^^^GTPyS or Cu^^^GTP chelation com- 
plex, and that Cu^^^nucleobase interactions further 
enforce Cu^^ selectivity and position the metal ion 
for Lewis acid catalysis. Replacing Mg^^ with 
[CoCNHaje]^^ significantly reduced product yield, 
but not kobsj indicating that the role of inner-sphere 
Mg^^ coordination is structural rather than catalytic. 
Replacing Mg^^ with alkaline earths of increasing 
ionic radii (Ca^^, Sr^^ and Ba^^) gave lower yields 
and approximately linear rates of product accumu- 
lation. Finally, we observe that reaction rates 
increased with pH in log-linear fashion with an 
apparent pKa = 8.0±0.1, indicating deprotonation in 
the rate-limiting step. 

INTRODUCTION 

Artificial ribozymes have been isolated by in vitro selection 
to catalyse a wide range of chemical reactions, although 



the chemical strategies used by the great majority of them 
are poorly understood. Revealing how they mediate these 
transformations is crucial for assessing the feasibihty of an 
RNA-based metabolism, for evaluating RNA world 
theories of early evolution, for engineering artificial 
enzymes and other tools for synthetic biology, and for bio- 
medical appHcations of ribozymes. Those ribozymes that 
have been studied in detail use several of the same catalytic 
strategies as protein enzymes, such as Lewis acid catalysis, 
proton transfer, precise positioning and desolvation of 
substrates and allosteric regulation, and a small handful 
of natural and artificial nucleic acid catalysts exploit the 
chemical reactivity of bound organic cofactors (1-4). 
Metal ion cofactors accelerate ribozyme catalysis by 
polarizing and acidifying inner sphere water, by increasing 
electrophilicity of phosphates and carbonyls, by stabilizing 
negative charges that develop on transition states, inter- 
mediates and products, by assisting in folding and by 
other mechanisms [reviewed in (5)]. Although Mg^^ is 
the dominant bioavailable divalent cation, the transition 
metal ions in the first row of the periodic table provide an 
array of unique chemical capacities, such as variations in 
Lewis acidity, charge density, preferred bond lengths 
and coordination geometries. Some metal ions retain 
their inner hydration sphere and interact via bound 
water molecules. For others, one or more waters are 
replaced by direct contact with the ribozyme 
or substrate(s). The relative stabilities of inner-sphere 
complexes involving these ions tend to follow the Irving- 
WilHams series (6): Mn^+ < Fe^+ < Co^+ < Ni^+ < 
Cu^+>Zn^+. A different series is observed for the 
affinities of these same ions in complexes with nucleotide 
mono-, di- and tri phosphates (7) and phosphate mono- 
esters (5), for which the affinities follow the order 
Mn^+ > Fe^+ > Co^+ ^ Ni^+ < < Cu^+ > > Zn^+ < Cd^+. 
Interestingly, Fe^^ has recently emerged as a potential 
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prebiotic metal ion, based on its likely high concentrations 
in aqueous solution on the early Earth and on its positive 
impact on ribozyme folding and catalysis (8). 

Ribozymes have been identified that use most of these 
ions. For example, the hammerhead (HH) ribozyme is 
active in low concentrations of Mn^^, Co^^, Ni^^, Zn^^ 
and Cd^^ (9-14), and RNaseP functions with Zn^^ as the 
sole divalent cation (15). From in vitro selections, RNAs 
have been isolated that assemble Ni^^ or Pt^^ into a 
coordinate-covalent RNA-amino acid complex (16) and 
that interact with immobilized Ni^^ or Zn ^ on affinity 
matrices (16,17). Use of divalent copper is highly 
unusual among structured nucleic acids — in particular 
among structured RNA — although a few examples have 
been noted among catalytic DNAs. A DNAzyme has been 
described that inserts Cu^^ and other transition metal ions 
into a mesoporphyrin ring (18,19), and DNAzymes have 
been identified that have strict dependence on Cu^^ for 
their catalysis of DNA ligation (20), DNA cleavage 
(21,22) and DNA-capping activities (23). Several 
5^-self-phosphorylating DNAzymes have been shown to 
use Cu^^ either as one among many functional cations 
[Dk5, (24)] or as a strict requirement [Cul, Cu4 and Cu7 
(25)]. Divalent ions are not always required, however, and 
the hairpin (HP), HH and Varkud Satellite (VS) ribo- 
zymes remain functional in the absence of divalent 
metals at high concentrations of monovalent ions, 
relying exclusively on nucleotide functional groups for 
catalytic chemistry (26-30). 

Proton transfer is integral to the mechanisms of ribo- 
zymes such as the hepatitis delta virus (HDV), HP, HH, 
glmS, VS and others (31-37). In the HDV ribozyme, for 
example, a bound metal ion hydroxide is in position to 
abstract a proton from the 2^0H nucleophile. The active 
site cytosine C75 is 'histidine-like', in that its apparent 
pKa is perturbed approximately three pH units from its 
normal value of 4.2 to neutrahty (32), and its N3 is in 
position to stabilize the leaving group by donating a 
proton to the ribose 5^0H of Gl (32-34). For self-cleavage 
by the HP ribozyme, the apparent pKa is ~6.0, and the 
deprotonated state of residue A3 8 has been shown to be 
important in catalysis (27,38,39). Intriguingly, in molecu- 
lar dynamics simulations, the Nl of A3 8 moves into close 
proximity of the active site 2^0H, where it would be in 
position to act as general base (40). For tertiary-stabilized 
HH ribozymes. X-ray crystallography show both G8 and 
G12 near the scissile bond, with G12 in position to act as 
general base (36,41), consistent with earher predictions 
based on pH rate profiles for site-specifically substituted 
ribozymes (35,42). The potential role of G12 as general 
base is further supported by labelling of its Nl via nucleo- 
philic attack on 2^-bromoacetamide (43). Nucleobase- 
metal ion interactions also appear to contribute to 
shifting nucleobase pKa values within the active sites of 
some extended HH ribozymes (9,44), and a suitably pos- 
itioned Mn^^ ion has been observed in HH ribozyme 
crystals (45). For the glmS ribozyme, an active-site 
guanine is in position to extract a proton from the 2^0H 
nucleophile (46-48), although it appears from studies of 
pH dependence of the reaction to participate in catalysis in 
its neutral, protonated state (49), perhaps by donating a 



hydrogen bond to the developing transition state. These 
observations for natural ribozymes contrast with studies 
of two previously studied kinase ribozymes (Kin. 46 and 
2PT3.2min) (50,51) and by a kinase DNAzyme (Dkl) 
(25). For all three catalysts, (thio)phosphoryl transfer 
rates were shown to be independent of pH over the 
range pH 6.0-8.0 (Dkl lost activity above 8.0). 

The present work details the metal ion and pH require- 
ments for phosphoryl transfer by ribozyme K28(1-77)C. 
This 58 nt RNA is a variant of ribozyme K28, which was 
originally selected as a 126 nt species for thiophosphoryl 
transfer activity using GTPyS as donor (52). Ribozyme 
K28(1-77)C folds into a compact pseudoknot and trans- 
fers a phosphoryl group from GTP (or a thiophosphoryl 
group from GTPyS) onto itself at two different sites in the 
primary sequence (Figure lA) (53). The original selection 
and functional analysis of ribozyme K28 and its deriva- 
tives were performed in the presence of alkaline earth 
metal ions (Mg^^ and Ca^^), transition metal ions 
(Mn^^ and Cu^^) and monovalent cations (K^ and Na^) 
and was buffered to near neutrality with hydroxyethyl- 
piperazine ethane sulphonate (HEPES, pH 7.5) (52). In 
the present work, we sought to determine the contribu- 
tions of each of these components to RNA-catalyzed 
phosphoryl transfer. We find that both outer-sphere and 
inner-sphere interactions with hydrated Mg^^ play im- 
portant structural roles. More surprisingly, K28(1-77)C 
is also the first kinase ribozyme to be fully dependent on 
Cu^^ for optimal activity, and the first for which the rate 
exhibits a log-linear pH dependence. Other ribozymes 
from the selection that gave rise to K28 and its derivative 
K28(1-77)C did not show this dependence on Cu^^. In-line 
probing (ILP) and competition studies with adenosine tri- 
phosphate (ATP) and ATPyS show no evidence of a direct 
Cu^^-RNA complex and instead reveal that Cu^^ enters 
into the ribozyme as part of a Cu^^»GTP chelation 
complex, with additional stabilization likely to come 
from interactions with nucleobase nitrogens in the RNA. 
The proposed binding mode explains the Cu^^ dependence 
and positions the bound Cu^^ to participate in the cata- 
lytic step of the reaction by serving as a Lewis acid 
catalyst. 

MATERIALS AND METHODS 

Materials 

Oligodeoxynucleotides were purchased from Integrated 
DNA Technologies (Coral ville, lA). RNA was transcribed 
in vitro using phage T7 RNA polymerase, which was 
overproduced in bacteria and purified in the laboratory. 
GTPyS and ATPyS were purchased from Sigma (St. 
Louis). Radiolabelled nucleotides for labelling internal 
phosphates or internal 2^0H positions ([oc^^P]-CTP and 
[y^^P]-GTP, respectively) were purchased from Perkin- 
Elmer (Waltham, MA). N-acryloyl-aminophenylmercuric 
(APM) chloride was prepared as described (53-55). 

Self-thiokinase reactions 

Ribozyme kinetic analyses were carried out essentially 
as described (53), in most cases using tri-layer 
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Figure 1. Cu^^ dependence of SP. (A) Secondary structure of 
K28(1-77)C. Circled sP indicate the two self(thio)phosphorylation 
sites; nucleotide Gl represents the major site (52). (B) Thiophosphoryl 
transfer from GTPyS was monitored while titrating Cu^^ from 0.01 to 
2[iM (0, 0.01, 0.1, 0.5, 1, 2, 5, 10, 50, 100 ^iM; fuU data set to 100 )iM is 
shown in inset.) The 'zero-Cu^^' data point was omitted from the plot 
to aUow the log plot. Maximal rate at saturation was 0.00151 ± 
0.00019 min~^ Fitting the initial rates to the quadratic form of a 
standard 2-state binding equation for a 1:1 complex yielded an 
apparent dissociation constant for Cu^^ (Kd^"^^) of 0.9±0.4)iM. 
(C) SP reactions (no sulphur) were performed with [y-^^P]GTP as phos- 
phoryl donor in SP buffer in the absence or presence of lOjiM Cu^^. 
Reaction times in hours are given above the lanes. 



organomercurial gels to analyse product formation by in- 
ternally radiolabelled transcripts (52,53,56). In brief, in- 
ternally radiolabelled K28(1-77)C RNA and other RNA 
molecules were gel purified, unfolded in water at 85°C for 
5min, then refolded on ice for 5min by addition of the 
same 5x self-phosphorylation (SP) buffer used originally 
to select this ribozyme in vitro (52) (Ix SP buffer = 6mM 
MgCl2, 0.2mM CaCl2, 0.5mM MnCl2, 10|iM CuCl2, 
200 mM KCl, 15mM NaCl, 25 mM HEPES, pH 7.5) or 
modifications of this buffer as detailed in the text. Kinase 
reactions used 1 |iM RNA (50 000-200 000 cpm) and were 
initiated by adding GTPyS to a final concentration of 
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0.5 to 1 mM and moving the reaction mixtures to 32° C. 
Except where noted, reactions were quenched at various 
times in stop buffer (95% formamide, 15mM ethylenedia- 
minetetraacetic acid and trace amounts of xylene cyanol 
and bromophenol blue as tracking dyes). 'Zero' time 
points were collected immediately after all components 
were present in the reaction mix. Products were separated 
on 8% denaturing tri-layered organomercurial gels in 
which the middle layer contained 100|ig/mL APM 
[(N-acryloylamino)phenyl]mercuric chloride (52,53,56). 
Autoradiographs were obtained with a FLA- 5000 
phosphorimager (FujiFilm) and analysed with 
MultiGauge software. The fraction of the RNA converted 
to product at a given time [f(t)] was calculated by dividing 
the intensities of RNA retained at the APM interface into 
the sum of all bands within a given lane. First-order ex- 
ponential rate constants (kobs) and extrapolated plateaus 
values (fmax) of most reactions were obtained by fitting the 
data to a first-order rate equation using KaleidaGraph 
(Synergy Software): /(O = (/max)[l - exp^-^^^^^>]. Data 
points in figure 5 are the mean values of at least two rep- 
licas for concentrations in which activity was detectable. 
Data points for the extrapolation of kobs in Figures 1 and 
5 were the mean values of at least three repHcas for con- 
centrations in which activity was detectable and were fit to 
a linear equation. 

Self-^^P-kinase reactions for Cu^^/sulphur dependence 

Ribozyme K28(1-77)C was transcribed without radiolabel, 
gel purified, unfolded in water at 85°C for 5min, then 
refolded on ice for 5min by addition of 5x SP buffer or 
5x SP buffer lacking Cu^^. Kinase reactions used 1 |iM 
RNA and were initiated by adding [y-^^P]GTP to a final 
concentration of 0.6 |iM and moving the reactions mixtures 
to 32°C. After the indicated times (0, 0.5, 1, 2, 3, and 9h), 
samples were moved to ice and ethanol precipitated imme- 
diately by addition of NaOAc, glycogen and ethanol to 
remove excess non-incorporated radiolabelled GTP. After 
resuspension in Ix gel loading buffer, samples were 
analysed on 8%, 8M urea denaturing polyacrylamide gels 
as aforementioned. 

Self-thiokinase reactions with competitors 

K28(1-77)C was unfolded and refolded as aforemen- 
tioned. ATP or ATPyS was added to a final concentration 
of 2mM and kept in ice for 5min. Donor GTPyS was 
added to final concentration of 0.5 mM, followed where 
indicated by addition of supplemental Cu^^ to a final con- 
centration of 30 |iM. 

ILP analysis 

Ribozyme K28(1-77)C was 5^-end labelled using 
[y-^¥]ATP, unfolded at 85°C for 5min, and then 
refolded by addition of SP buffer with or without CuCl2, 
and modified to contain 25 mM TRIS (2-Amino- 
2-hydroxymethyl-propane-l,3-diol) (pH 8.0) in place of 
25 mM HEPES (pH 7.5). Competitors, donor and supple- 
mental Cu^^ were added in the same order as described 
earlier in the text. ILP reactions were carried out at room 
temperature and stopped after lOh by addition of half 
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a volume of 95% formamide and 50 mM ethylene- 
diaminetetraacetic acid. Digestion products were then 
separated on an 8% denaturing polyacrylamide gel electro- 
phoresis gel, and autoradiographs were analysed as 
aforementioned. 

Self-thiokinase reactions for pH dependence 

SP buffer was modified as follows: 25 mM [2-(A^- 
morphoHno)ethanesulphonic acid] (MES) was used in 
place of HEPES for pH 5.5, 6.0. 6.5, 6.8; HEPES was 
used for pH 6.8, 7.0, 7.2, 7.5, 7.8, 8.0; and TRIS was 
used at pH 8.0, 8.2, 8.5, 8.8, 9.0, 9.5. Owing to precipita- 
tion of manganese at elevated pH, this cation was 
eliminated in all SP buffers used in the pH study. 
Control reactions performed over the pH range 5.5-8.0 
in HEPES verified that activity remains the same with 
and without Mn^^ (data not shown), thereby validating 
this approach. Very low product accumulation at the 
lowest pH values (MES buffer) precluded fitting to a 
first-order exponential rate equation. Rates for these low 
pH values were obtained from the slope of the Hne of best 
fit through these data, adjusted to an assumed plateau of 
66% (plateau values for the HEPES and TRIS data are 
66 ± 7% conversion to product). The observed rate con- 
stants were then fit to a standard equation for a 
one-proton transfer: kobs = kmax/[l + lO^^^^^'^^^]. Kinetic 
reactions were performed at least in duplicate. 

RESULTS 

A Cu^^-dependent kinase ribozyme 

Ribozyme K28, which is the 126-nt parent form of the 
58-nt ribozyme K28(1-77)C, was originally selected for 
thiophosphoryl transfer activity in SP buffer, which 
contains monovalent ions (200 mM and 15mM 

Na^), divalent alkaline earth ions (6mM Mg^+ and 
0.2 mM Ca^^) and divalent transition metal ions 
(0.5 mM Mn^^ and 10 |iM Cu^^). To begin to understand 
the contributions of each of these ions to thiophosphoryl 
transfer activity by ribozyme K28(1-77)C, reactions were 
carried out in SP buffer in which individual components 
were omitted. Surprisingly, the ribozyme was essentially 
inactive without Cu^^. When product formation was 
monitored, as Cu^^ was titrated from 0.01 to 100 |iM in 
the presence of all of the other buffer components, initial 
rates increased quickly at sub-micromolar concentrations 
and reached a plateau near the concentration of the RNA 
(1 |iM), with Httle additional change at still higher concen- 
trations (Figure IB). Fitting the initial rates to the quad- 
ratic form of a standard two-state binding equation 
yielded a good fit for a 1:1 stoichiometry and an 
apparent dissociation constant for Cu^^ (Kd^^^^) of 
0.9 ± 0.4 |iM, which was close to the concentration of 
RNA used in the assay. 

Cu^^ is considered to be 'borderline' on the Pearson 
scale of hard and soft metal ions, and it associates 
well with 'soft' ligands such as nitrogen and sulphur (57). 
To determine whether the observed Cu^^ requirement 
is the result of essential interactions between the Cu^^ 
and the thiophosphoryl group of the GTPyS donor, 



non-radiolabelled K28(1-77)C RNA was incubated with 
[y-^¥]GTP in SP buffer with and without Cu^+. 
Accumulation of radioactivity in the RNA under these con- 
ditions indicates sulphur-independent SP by K28(1-77)C. 
As with the GTPyS donor, product formation from the 
OTP donor occurred much more vigorously when Cu^^ 
was present, even though there is no sulphur in the GTP 
donor (Figure IC). Therefore, the observed requirement for 
Cu^^ in both GTP-dependent and GTPyS-dependent reac- 
tions indicates a specific function of the metal ion and not a 
spurious consequence of having used GTPyS during the 
initial selection. Therefore, all further reactions included 
1 0 |iM Cu^^ unless otherwise noted. 

Uniqueness of the Cu^^ dependence of kinase 
ribozyme K28 

We next sought to determine whether other divalent ions 
could substitute for Cu^^ and whether other ribozymes 
from the same selection as K28 displayed a similar de- 
pendence on Cu^^. Little or no product was observed in 
overnight reactions when 10|iM Cu^^ was replaced with 
this same concentration of transition metal ions (Cr^^, 
Co^+, Ni^+, Cd^+, Mn^+ and Zn^+) or alkahne earth 
metals (Mg^^, Ca^^, Sr^^ and Ba^^), in the presence of 
the other buffer components (Figure 2, data not shown). 
These other metal ions are therefore poor substitutes for 
10|iM Cu^^. Because Cu^^ binds phosphate more tightly 
than do these other metal ions (5,58), we next monitored 
product yields as a function of metal ion concentration. 
There is a modest increase in product formation when 
10 uM Cu^+ is replaced with at 100-500 |iM Mn^+ or 
Ni ^, but strong inhibition is observed for Zn^^, Co^^ 
and potentially Cd^^ at >100|iM concentrations, 
perhaps because of non-specific binding to the RNA. 
Therefore, the role played by Cu^^ cannot be fully 
provided by these other divalent metal ions, which can 
only support sub-optimal activity (at best) even at 10- to 
100-fold higher concentrations. 

The in vitro selection that gave rise to kinase ribozyme 
K28 also produced several other ribozymes that could be 
classified into at least three different families based on 
their secondary structures (52). Because all of these ribo- 
zymes were selected under identical ionic conditions, two 
ribozymes from each structural family were tested for 
Cu^^ dependence by performing the self-thiphosphory- 
lation reaction in the absence and presence of 10|iM 
Cu^+. Although ribozymes K28(1-77)C and K28 showed 
the expected dependence on Cu^^, all six of the other ribo- 
zymes produced at least as much product in the Cu^^- 
depleted SP buffer as they did in the complete SP buffer, 
and K5 and K6 yielded slightly more product when Cu^^ 
was omitted (Figure 3). Cu ^ dependence is therefore not a 
widespread characteristic among the ribozymes that were 
selected in the presence of Cu^^ and is unique for K28 and 
its derivatives. 

ILP identifies RNA binding by GTP and GTP-Cu^^, but 
not by Cu^^ alone 

ILP assays were used to determine whether ribozyme 
K28(1-77)C uses Cu^^ to fold into its active structure. 
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Figure 2. Substitutions for Cu ^ reduce or eliminate activity. Product yields were plotted for 18 h self-kinase reactions that were carried out in SP 
buffers that were modified by replacing 10|iM CuCl2 with MnCl2, NiCl2, ZnCh, CdCl2 or C0CI2 at the indicated concentrations (grey bars). Yield 
from a parallel reaction carried out in normal SP buffer ('lOjiM Cu^^') is plotted for comparison (black). Average values for four sets of reaction 
(except for Cd^^, for which N = 2) are plotted. Error bars represent the standard error of the mean. 



ILP gives information on conformationally mobile 
phosphoester bonds that increase the fraction of time in 
which the 2^0H samples conformations appropriate for 
in-line attack on the adjacent phosphate and release of 
the downstream 5^0H (59). No new ILP cleavage was 
observed when Cu^^ was added to the reaction buffer at 
a concentration of either 10 or 30 |iM (Figure 4A, lanes 
1-3). In contrast, addition of 0.5 mM GTP induced a 
strong cleavage after A32, and this cleavage was further 
sensitized by addition of Cu^^ (Figure 4 A lanes 4-6). No 
changes are seen at other positions in response to addition 
of Cu^^ or GTP. These data support direct binding of 
GTP in the absence of Cu^^ and enhanced binding of 
GTP when Cu^^ is present, but they do not provide 
evidence of a direct interaction between RNA and Cu^^ 
in the absence of GTP. Instead, we interpret the stimula- 
tion as indicating that Cu^^ binds the triphosphate of 
GTP in solution, and that the metal-bound form of the 
GTP makes additional RNA-Cu^^ contacts. 

Cu^^ interacts with the triphosphate region of the 
phosphoryl donor 

Previous studies measured competition for access to the 
ribozyme active site between GTPyS and a series of 



substrate analogs (52). GTP, 3^dGTP and a few other ana- 
logues were recognized by the RNA and occupied the 
active site in a manner that prevented GTPyS use. ATP 
exhibited no such competition and is therefore not 
beheved to interact with the ribozyme. Consistent with 
this interpretation, the intensity of ILP cleavage at A32 
in the presence of 2mM ATP and 0.5 mM GTP was iden- 
tical to those observed in reactions without ATP, and 
addition of excess Cu^^ had no effect on the cleavage in- 
tensity (Figure 4 A, lanes 7 and 8). However, we reasoned 
that a different mode of competition could result from 
sequestration of free Cu^^ by the triphosphates of other 
nucleotides, thereby reducing both GTP-induced ILP 
cleavage and self-thiophosphorylation yield from the 
GTPyS donor. Furthermore, although Cu^^ dependence 
does not require the sulphur moiety of GTPyS (Figure 2), 
the presence of the sulphur can be safely assumed to 
mcrease affinity of the NTP for Cu^+. Indeed, when ILP 
was carried out in the presence of 2mM ATPyS, the 
cleavage signal at A32 owing to GTP was diminished. 
Signal strength was partially restored on supplementation 
with additional Cu^^ (Figure 4A, lanes 9 and 10). 
Quantification of signal strength at position A32 
shows that ATPyS reduced total cleavage by ~4. 6-fold 
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Figure 3. Stimulation by Cu^^ is unique to ribozyme K28. Percentage 
of the indicated ribozymes that were converted to 
self-thiophosphorylated product is shown for 18-h reaction that was 
performed in the presence (black) or absence (grey) of lOjiM Cu^^. 
The secondary structural features shared between each pair of ribo- 
zymes are indicated above the graph [detailed in ref. (52)]. Asterisks 
indicate sites of self-thiophosphorylation; '28C' indicates ribozyme 
K28(1-77)C. 



(Figure 4B). When GTP was omitted from the ILP 
reaction, neither ATP nor ATPyS induced RNA cleavage 
at A32, irrespective of Cu^^ concentration (Figure 4A, 
lanes 11-14). These results rule out the possibility of 
direct RNA-ATP or RNA-ATPyS interactions and 
suggest that competition by ATPyS is indirect, as would 
be expected from sequestration of free Cu^^ by ATPyS. In 
the previous analogue studies, GMP competed only weakly 
with GTPyS, indicating that the beta and gamma phos- 
phates are important for efficient recognition. ILP reac- 
tions using GMP in place of GTP gave a Cu^^-dependent 
cleavage signal at A32 that was much weaker than the 
GTP-dependent signal, along with faint bands at C39, 
A46 and C47, all of which were suppressed by ATPyS 
(Supplementary Figure SI). 

The effect of ATPyS on the ribozyme' s abihty to 
interact with GTPyS was also examined by monitoring 
self-thiophosphorylation. Ribozyme K28(1-77)C was 
incubated overnight with 0.5 mM GTPyS and 2.0 mM 
ATP or ATPyS. As observed previously (52), excess 
ATP did not compete and had no effect on reaction 
yield (Figure 4C). In contrast, ATPyS reduced the 
observed self-thiophophorylated product by ~3. 8-fold, 
and the addition of supplemental Cu^^ restored the yield 
of thiophosphorylated product. These data estabhsh that 
Cu^^ does interact with the triphosphate region of the 
nucleotide donor, and that the effects of ATPyS on ILP 
cleavage and thiophosphoryl transfer is due to chelation 
of Cu^^ by the thiophosphate of ATPyS. 

Non-specific Mg^^ interactions 

Apart from Cu^^, the major divalent ion present during 
the initial selection was Mg^^. Therefore, all divalent ions 
except for Cu^^ were omitted from SP buffer, and 
self-thiophosphorylation was monitored for ribozyme 
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Figure 4. Triphosphate interactions revealed by ILP. ILP analysis was 
used to determine Cu^^-induced structural changes. (A) ILP reactions 
were performed for 10 h in SP buffer modified to contain 25 mM TRIS 
(pH 8.0). ATP or ATPyS was added before GTP when both nucleotides 
were present. For lanes 3, 6, 8 and 10, Cu^^ concentration was initially 
lOjiM and was adjusted to 30)iM after addition of ATP/ATPyS and 
GTP. (B) Integrated ILP cleavage signal at position A32 under the 
indicated conditions. (C) Chelation-mediated competition in self-thiopho- 
sphorylation reactions was evaluated by quantifying product formed after 
18 h reactions under the indicated conditions. As aforementioned, com- 
petitors were added before addition of GTP, and supplemented Cu^^ was 
added after the addition of competitors and donor. 



K28(1-77)C, as Mg ^ concentration was titrated from 
0.01 to 15mM. Thiophosphoryl transfer rates increased 
linearly from 1 to 15mM Mg^^, with no indication of 
approaching saturation over this range (Figure 5A). 
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Figure 5. Outer-sphere and inner-sphere Mg^^ interactions for ribozyme K28(1-77)C. (A) Reactions were carried out in the presence of various Mg^^ 
concentrations (0, 0.01, 0.05, 0.1, 0.2, 0.5, 1, 3, 6, and 15 mM). Observed kinetic constants (kobs) are plotted versus Mg^^ concentration for reactions 
in which product was detected. Reactions in which Mg^^ was < 1 mM yielded insufficient product for accurate determination of kobs and were not 
plotted. (B) Accumulation of thiophosphorylated product in the presence of the indicated Mg^^ concentrations after short (3h) reactions in the 
presence (black) or absence (grey) of lOjiM Cu^^. (C) To evaluate potential inner-sphere Mg^^ interactions, kinetic traces were plotted for reactions 
in which the 6mM Mg^^ component of SP buffer was replaced with 6 or 12mM CoHex. Uncertainty reported for the rate constants reflects the 
uncertainties of the fit. (D) Similar to part C, except that Mg^^ was replaced with other alkaline earth metal ions. Dashed fines indicate linear fits to 
these data. 



Titration with Mn yielded a similar pattern as the Mg 
titration to ~6 mM but was inhibitory at still higher con- 
centrations (data not shown). Little or no product was 
observed with < 1 mM Mg^^ or Mn^^, and no product 
was detected when 6mM Co^^ or Ni^^ was used in 
place of 6mM Mg^^ (data not shown). Thus, Mg^^ or 
Mn^^ (or potentially other divalent ions) is required for 
phosphoryl transfer by K28(1-77)C via low-affinity inter- 
actions, but neither divalent cation is specifically required. 

Natural kinase and polymerase enzymes make extensive 
use of Mg^^^NTP and Mg^^»dNTP complexes [see 
(60-62) and references therein]. However, Cu^^ has a 
higher affinity for phosphates than do Mg^^ and other 
divalent metal ions [> 30-fold higher for NMP-metal ion 
complexes (5)]. Formation of Cu^^»GTPyS is therefore 
expected to be strongly favoured over formation of 
Mg^^^GTPyS when both ions are present. To test 
whether Cu ^ dependency at relatively low Mg^^ concen- 
trations could be overcome at higher Mg^^ concentra- 
tions, self-thiophosphorylation yield for relatively short 
reaction times (3h) was monitored at Mg^^ as high as 
lOOmM in both the presence and absence of 10|iM 
Cu^^. As aforementioned, product yield in the presence 
of Cu^^ continued to increase with Mg^^ concentration 
with no sign of approaching saturation (Figure 5B, 



black bars). Interestingly, product yield also increased to 
well above background in the absence of Cu^^ at elevated 
Mg^^, although the yield reached a plateau at ~30mM 
Mg^^ and did not approach the yield observed in the 
presence of Cu^^ (Figure 5B, grey bars). Thus, although 
Mg^^ may partially substitute for Cu^^ at high concentra- 
tions (for example, by forming a Mg^^^GTPyS complex), 
this partial rescue is only observed when Mg^^ is least 
3000 times higher than the Cu^^ concentration and >15 
times higher than the normal Mg^^ concentration in SP 
buffer. It would therefore be premature to infer that the 
precise contribution of Mg^^ at high concentrations is 
fully interchangeable with the role of 10|iM Cu^^. 
Overall, the observations from Mg^^ titrations suggest 
that low-affinity, non-specifically bound Mg^^ ions serve 
primarily to aid folding the ribozyme into its active struc- 
ture, which can then use the Cu^^^GTPyS complex for 
optimal catalysis. 

Evidence for inner sphere interactions with Mg^^ 

Mg^^ can contribute to ribozyme activity by contacting 
ribozymes and substrates directly (inner sphere), and via 
ligands such as water that are bound directly to the metal 
ion (outer sphere). Cobalt hexammine {[Co(NH3)6]^^, 
CoHex} is a reasonably good analogue for outer sphere 
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interactions with hydrated Mg owing to its similar size, 
geometry and charge, but kinetic stabiHty of the coord- 
inated NH3 Hgands precludes significant inner-sphere 
interactions. CoHex supports catalysis as well as, or 
nearly as well as, Mg ^ for ribozymes that do not 
require direct contact with the metal ion, such as the HP 
ribozyme (26,27,63) and the in vitro selected ATRib/AT02 
acyltransferase ribozymes (64,65). In contrast, CoHex 
does not support catalysis at all, or not nearly as well as 
Mg^^, for ribozymes that exploit direct contacts with the 
metal ion, such as the HH (44,66), VS (67), RNaseP (68), 
HDV (69,70) ribozymes and the in vitro selected 2PT3.2 
kinase ribozyme (51), and CoHex even induces aberrant 
folding of RNaseP (15). 

When the 6mM Mg^^ component of SP was replaced 
with 6mM CoHex, self-thiophosphoryl transfer was 
greatly reduced (~5% yield after 16 h) (Figure 5C) 
relative to Mg^^ at the same concentration, supporting a 
potential role for at least one inner-sphere contact with 
Mg^^. DoubHng the CoHex concentration increased 
yield after 16 h (18%), but the calculated kobs values at 
these two CoHex concentrations are similar (0.0034 ± 
0.0006 and 0.0026 ± 0.0003 min"\ respectively). Increas- 
ing CoHex therefore stimulates the amount of RNA that 
can self-thiophosphorylate without altering the intrinsic 
rate. Importantly, the kobs values measured in CoHex/ 
Cu^^ are not significantly different from those measured 
in Mg^+/Cu^+ (0.0032 ± 0.0001 min"^), even though much 
more product is formed in Mg^^/Cu^^. Thus, Mg^^ plays 
a role in forming the active structure of ribozyme 
K28(1-77)C but does not play a catalytic role in 
thiophosphoryl transfer, even for the Mg^^ ions that 
appear to make inner sphere contacts. 

Probing the dimensions of the Mg^^-binding site(s) 

To determine whether divalent alkaline earth ions larger 
than Mg^^ could also support catalysis, thiophosphoryl 
transfer was monitored for reactions in which Ca^^, Sr ^ 
or Ba^^ was the only divalent ion apart from 10 |iM Cu^^. 
Replacing 6 mM Mg^^ with 6 mM Ca^^ or Sr^^ produced 
a modest reduction in overall rate and yield during the 11- 
h reaction, and much less product formed in 6 mM Ba^^ as 
the major divalent cation (Figure 5D). Interestingly, the 
rate of product accumulation for these ions is approxi- 
mately linear rather than fitting well to a first-order 
kinetic best-fit curve, potentially indicating a shift in the 
rate-limiting step (e.g. slow folding). These three trends 
(reductions in rate and yield, and deviation from simple 
first-order kinetics) all correlate with the ionic radius of 
these alkahne earth metal ions. Specifically, thiopho- 
sphoryl transfer activity decreased in the order 
Mg^^>Ca^^>Sr^+»Ba^+, for which the ionic radii of 
their hexacoordinate forms are 0.86 A, 1.14 A, 1.32 A 
and 1.49 A, respectively (71). 

pH-dependence 

To determine the potential role of proton transfer in the 
reactivity of ribozyme K28(1-77)C, apparent rate con- 
stants were measured, as pH was varied from pH 5.5 to 
9.5. Three different buffers were used to span this range 



(MES, HEPES and TRIS), including separate reactions in 
MES or HEPES (pH 6.8) and in HEPES or TRIS (pH 8.0) 
at the two overlapping pH values to control for 
buffer- specific effects. In contrast to previously analysed 
kinase ribozymes (50,51) and DNAzymes (25), for which 
the rates showed no pH dependence, thiophosphoryl 
transfer by K28(1-77)C was strongly dependent on pH 
(Figure 6). The observed first-order rate constants (kobs) 
increased ~ 1000-fold from pH 5.5 to 8.5, indicating one 
proton-transfer event in the rate-limiting step. Fitting 
these data to a single-proton equation (see 'Materials 
and Methods' section) yielded an apparent pKa value of 
7.99 ± 0.08, and a maximal rate constant 
(0.049 ± 0.002 min"^) that is > 15-fold higher than that 
observed in normal SP buffer. The rate constants for re- 
actions carried out in HEPES were 3-5-fold lower than 
those for reactions carried out in either MES or TRIS at 
the two overlapping pH values, and the overall increase 
with pH was less steep in the presence of HEPES relative 
to MES or TRIS, indicating buffer-specific effects that 
make HEPES a suboptimal buffer for ribozyme 
K28(1-77)C. 



DISCUSSION 

We demonstrate here that phosphoryl transfer by the 
K28(1-77)C kinase ribozyme strongly depends on the 
presence of Cu^^, that the reaction also requires Mg^^ 
or similar divalent cations and that the rate-limiting step 
of the reaction includes a proton transfer event. Very Httle 
product was formed in the absence of Cu^^, independent 
of whether the donor was GTPyS or [y-^¥]GTP, estab- 
lishing that the Cu^^ requirement persists apart from any 
potential interaction with the donor sulphur. Fitting the 
titration data to a binding isotherm indicated a 1:1 stoi- 
chiometry, with Kd^''^^ of 0.9 ± 0.4 |iM. However, ILP 
analysis argues against interpreting these data in terms 
of a simple direct metal ion-RNA interaction, as there 
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Figure 6. pH rate profile of K28(1-77)C. kobs values are plotted as a 
function of pH. Because HEPES (triangles) exhibited sub-optimal, 
buffer-specific effects, only the data from reaction carried out in 
MES (circles) and TRIS (squares) were used in determining the fit. 
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Figure 7. Schematic representation of metal ion interactions in 
ribozyme K28(1-77)C. Cu^^ (blue) is proposed to serve as Lewis acid 
catalyst by chelating both the nucleotide triphosphate and nucleobase 
nitrogens within an azophilic binding pocket in the RNA, while Mg^^ 
(red-orange) contributes to structure through both inner sphere and 
outer sphere interactions. Postulated RNA-GTP interactions through 
the guanine Watson-Crick face, 3'OH and phosphates are based on 
previous analysis using analog competition (52). 



are no new RNA cleavages associated with addition of 
Cu^^ alone. Instead Cu^^ stimulates GTP-dependent 
RNA cleavage and partially alleviates the inhibitory 
effects of competitor ATPyS on RNA cleavage during 
ILP and on product formation during self-thiophosphory- 
lation reactions. All of these observations suggest that 
Cu^^ enters the ribozyme as a Cu^+#GTPyS chelation 
complex, wherein it is additionally chelated by a 
nitrogen-rich site in the RNA (Figure 7). Cu^^ is highly 
azaphihc and prone to Jahn-Teller geometry distortions in 
certain binding environments. As such, it may be particu- 
larly well suited to bind within a nitrogenous pocket that 
is unable to form strong interactions with other metal 
ions. Such interactions would ideally position the metal 
ion to promote the chemical step of the reaction by with- 
drawing electrons from the phosphate (Lewis acid cataly- 
sis), making it more susceptible to nucleophilic attack by 
the acceptor oxygen. In addition, Cu^^ binds phosphate 
more tightly than do the other metal ions assayed here 
(5,58), and phosphodiester cleavage rates by HH ribo- 
zymes in various metal ions correlate with binding 
affinities of those metal ions with phosphate monoesters 
(5). However, phosphate affinity alone does not explain 
metal ion selectivity, as these other metals cannot fully 
substitute for Cu^^ even at 100-fold (for transition metal 
ions) to > 1000-fold (for Mg^^) higher concentrations. 
Partial rescue in the absence of Cu^^ at very high Mg^^ 
concentrations suggests that the Mg^^^GTPyS complex 
can also be used. However, the yield of thiophosphory- 
lated product in the presence of Cu^^ continues to increase 
well above the Mg^^ concentration where the yield satur- 
ates in the absence of Cu^^, indicating that the 
Mg^^^GTPyS is not used by the ribozyme as well as the 
Cu^^»GTPyS complex. We therefore propose that Cu^^ 
ion selectivity results from a combination of coordination 
geometry, azaphilicity, geometric constraints of the active 
site and affinity for phosphates. 

Apart from K28(1-77)C, there are few or no other Cu^^- 
dependent ribozymes available for comparison. A Diels- 
Alderase ribozyme selected in the presence of Mn^^, Fe^^, 



Co^^, Ni^^, Cu^^, Cd^^ and Zn^^ was initially reported as 
requiring Cu^^ (72,73). They interpreted the Cu^^ require- 
ment in terms of catalytic Cu^^ Lewis acid sites, consistent 
with Lewis acid-catalysed Diels-Alder reactions in water. 
However, one of those Diels-Alderase ribozymes was later 
shown to function in the absence of Cu^^, given elevated 
concentrations of Mg^^, thereby estabhshing that Cu^^ 
plays a non-essential role in that ribozyme (74). 
Nevertheless, Lewis acid catalysis figures prominently in 
ribozyme reactions related to phosphoryl transfer. For 
example, a recent crystal structure of the HDV ribozyme 
shows inner-sphere ligand interactions between an active 
site Mg^^ ion and the pro-R(P) oxygen of the scissile phos- 
phate and the 2^-hydroxyl nucleophile (75). Similarly, 
polymerase and pyrophosphatase enzymes (which 
formally catalyse phosphoryl transfer to a polynucleotide 
chain or to water) use two Mg^^ ions as Lewis acid cata- 
lysts to shield charges and lower the pKa of the leaving 
group (62). 

The reduction in overall thiophosphoryl transfer 
activity by ribozyme K28(1-77)C on replacing Mg^^ with 
the non-exchanging CoHex indicates that the ribozyme 
takes advantage of inner sphere interactions with Mg^^; 
however, the primary role of Mg^^ is to help estabUsh the 
active structure, as evidenced by the insensitivity of the 
first-order rate constants to the CoHex substitution. 
Larger alkaline earth metal ions also promoted 
self-thiophosphorylation, but they did so less well and 
induced an intriguing and increasingly prominent devi- 
ation from first-order kinetic behaviour. The correlation 
of these trends with ionic radius may reflect simple steric 
constraints (upper size limit on ions that can successfully 
replace Mg^^ in at least one specific binding site), or they 
may arise from chemical features of the ions that also 
correlate with ionic radius, such as increased polarizabiHty 
or reduced charge density. 

In addition to its unusual metal ion requirements, 
K28(1-77)C is also the first kinase ribozyme to demon- 
strate a dependence on pH. Prior work with ribozymes 
Kin. 46 (50) and 2PT3.2min (51) observed no change in 
reaction rates, as pH were varied over several pH units. 
These observations suggested that deprotonation and ac- 
tivation of the 2^0H acceptor nucleophile was not rate 
limiting for those ribozymes and potentially suggested 
an SNl-like dissociative mechanism in which release of 
the gamma (thio)phosphate from the donor precedes for- 
mation of the acceptor-phosphate bond (50). The 
log-Hnear increase in the pH rate profiles observed here 
indicates that deprotonation of a species with pKa of 7.99 
± 0.08 is required for maximal thiophosphoryl transfer 
activity. The ribose 2^0H at the two modification sites are 
logical candidates for deprotonation to increase their 
nucleophilicity for attack on the gamma phosphate. 
However, the pKa for a ribose 2^0H is in the range of 
12.2-13.7 (52,76,77) and would need to be perturbed four 
to five pH units to account for the data here. The pKa for 
hydrated Mg^^ (11.4) is closer to the observed kinetic 
pKa, but the essentially unperturbed kobs in CoHex 
argues against a catalytic role for Mg^^. The pKa values 
for free guanosine (N1)H and uridine (N3)H are each 
~9.2 and could be shifted to near 8.0 within the context 
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of the ribozyme active site, for example, via metal ion 
interactions (9,44). Intriguingly, the pKa of Cu^^-bound 
water is near 7.5 (71,78) and is expected to shift upward in 
complexes with ligands that donate electrons to the metal 
ion centre. Deprotonation of one of the waters on a bound 
(inner sphere, partially dehydrated) catalytic copper 
hydrate or Cu^^-nucleobase interactions that modulate 
the pKa of the nucleobase would provide an intriguing 
functional linkage between the Cu^^ requirement and 
pH sensitivity of this kinase ribozyme. 

Although HEPES is generally considered to be a 
non-chelating buffer with respect to hard metal ions 
such as Mg^^ (79), interactions between Cu^^ and 
HEPES have been observed by isothermal titration calor- 
imetry (80), electron paramagnetic resonance spectros- 
copy (81) and potentiometric titrations (81,82). 
However, the suboptimal performance of HEPES buffer 
is not fully explained by a simple competition between the 
buffer and the ribozyme/donor for access to the metal ion; 
TRIS buffer forms TRIS#Cu^+ chelates with affinities that 
are similar to those of the HEPES»Cu^^ complex (58,81) 
and yet is well behaved with respect to its influence on the 
pH dependence of the reaction rate. Buffer-specific effects 
were not further explored. 

In conclusion, the chemical mechanism of the ribozyme 
K28(1-77)C stands in sharp contrast with previously 
characterized kinase ribozymes, using a tightly bound 
Cu^^ ion and also using a species that deprotonates with 
pKa ^ 8.0, in addition to specifically bound (potentially 
inner sphere) and non-specific cations required for folding. 
The contributions of Cu^^ to the catalytic mechanism 
include unique chemical features that cannot be 
provided by other first-row transition metal ions and 
that appear to implicate Lewis acid catalysis through 
Cu^^ chelation by the incoming donor triphosphate in a 
Cu^^^GTP or Cu^^^GTPyS complex. The chemical 
strategies that ribozymes K28 and its derivatives use to 
accomphsh (thio)phosphoryl transfer are clearly distinct 
from those used by other kinase ribozymes to perform 
the same reaction, including those that were co-selected 
along with K28 under identical ionic conditions. An im- 
plication for RNA world theories is that the parallel evo- 
lution of strikingly diverse metal ion and pH requirements 
among ribozymes within a single population highUghts the 
ability of RNA to evolve structures that capitalize on a 
diverse suite of mechanisms and local resources to accom- 
plish the same chemical objectives. 
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